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[дентификащя безперервно! функци 
в одновим!рному парабол!чному равнянн! 


В статье рассматривается задача идентификации непрерывной функции параболического уравнения в 
частных производных. Найдено аналитическое выражение для расчёта градиента неявно заданного 
функционала. Для определения градиента используется модернизированный классический метод 
множителей Лагранжа. 
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Розглядаеться задача 1дентификацй параметра у вигляд! безперервно! функцй парабол!чного рвняння в 
приватних похдних. Знайдено аналтичне вираження для розрахунку гращента неявно заданого функцо- 
нала. Для визначення градента використовуеться модернзований класичний метод множникв Лагранжа. 
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Пусть в пространственно-временной области > = [х,,х,]х[#,#,] функция Х(х,0) 


а?’Ь 
удовлетворяет квазилинейному параболическому уравнению: 
д д д 
5% я [6 9 9 - 0, (12), (1) 


Ох 


где «(/), В(х,г) непрерывные или кусочно-непрерывные функции. 
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Будем считать, что в области 1 =[х0,х1]х[1,,%ь] функция Д(х,1) = (1), а в области 
У =У\51 =[хр,х2]х[1.,] функция Д(х,!) задана и принимает значение р/. Поэтому 
можно записать 
Вх.) = В + [В (В 9(-ху), (2) 
где В, - непрерывная идентифицируемая функция. При этом уравнение (1) 
преобразуется к виду: 


2 
а Ру°-- Вод, ей =0, (1,х)ЕХ. (3) 
1 ох Ох 
Г раничные условия зададим в виде: 
о =0, А =-Иу-Я), (4) 
й х=х0 х х=хо 


гдеу, Я — некоторые константы. 
Начальное условие имеет вид: 
ах) = Е, (ьхе [ко,х2 |. (5) 
Качество идентификации функции Вь(®) оценим функционалом: 
1 
ИВь)= [хе - рахоРа, (6) 


‘а 


где Л, — экспериментально определенная функция /(1,х.), х. ®[хо»х|). 
Сформулируем задачу идентификации следующим образом. Необходимо найти не- 
прерывную функцию 2,(!) в уравнении (3), которая доставляет минимум целевому 


функционалу /(Дь). 

Наиболее эффективными методами идентификации являются прямые экстремаль- 
ные методы [1-4]. Они используют градиент целевого функционала для итерационных 
коррекций идентифицируемого параметра. 

Получим аналитическое выражение градиента целевого функционала ./(р,). 
Методом определения градиента в настоящей работе, как и в работах [4-6], является 
модернизированный классический метод множителей Лагранжа [2]. Данный метод вклю- 
чает в себя следующие этапы: 1) линеаризация уравнения (1) и целевого функционала (6) 
относительно д }еЕТ(»), д ВеИ(>); 2) отображение линеаризованных уравнений и функ- 
ционала в пространство К; 3) преобразование получившихся отображений вариаций к 
отображениям независимых вариаций 5 и 0рВ; 4) объединение элементов задачи иден- 
тификации в одинаковых пространствах; 5) выделение градиента целевого функционала. 
Через ГС) и ИС.) обозначены пространства состояний и управлений, соответственно, 


ГС), ИС) Е [> ‚ где Г — пространство функций с интегрируемым квадратом. Перейдем к 


определению градиента У.Л. 
Уравнение (3) для дальнейших преобразований, удобно записать в виде сис- 


темы с формальной переменной 4 = В(х,1 р : 
и аа 
[61 Ох 7 
5 (1) 
а4-В— =0, (.х)ЕХ. 
Ох 
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Система уравнений, линеаризованная относительно д }еЕТ(>), дВеИ(>), имеет 


вид: 
09 д 
ау 9-0 7 
я Ох га (8) 
ов а =0 ЕГ(>). 
Граничные и начальные условия имеют вид: 
М =0, 4 _ =хд = обл —79 , (9) 
9.) =0. (10) 
Линеаризованный функционал (6) принимает вид: 
1Ь 
&7= [2 Рб,х,)- ху = 2,9 ЕК, |5 
а 


где звездочка над символом обозначает пространство, сопряженное к исходному. 
Для отображения линеаризованной системы (8) в пространство К введем ли- 


нейный функционал-вектор © = =, 82 Е Р (>) . Умножим скалярно данный функцио- 


нал на вектор е = _ 
г 
1ъх2 
(ве) = || [(впел+ вое и =0ЕК. (12) 
'ах0 
Преобразуем выражение (12) к виду скалярного произведения относительно ва- 
риаций 6, да. Для этого необходимо преобразовать следующие слагаемые: 


16) д 
а г) 79- = 27 289. _ ау"; 
вв = ИРУ -дыл 5. 


Полученные дополнительные слагаемые в виде производных легко интегри- 
руются по х и {Е ввыражении (12). Окончательно, с учетом (9), (10), получаем: 


1х2 
(в.е)= | (-= ЕЕ а: ИЕ ОЕ м $ 


1а%0 
1 
к оовы, &- еооваа, &- |(с19 + Вво9Г) „+ 


хо хо ‘а 
1 
Е [ст + Азат) =0 ед. (13) 
12 
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Теперь можно объединить выражение (13) с линеаризованным функционалом (11). 
Для того чтобы избавиться от компоненты, принадлежащей сопряженному прост- 


ранству ис) ‚ потребуем: 


де д | 
о - 72 у 1.0 ву", 

ГА дх (14) 
а е* (5). 


Исключая из уравнений (14) линейный функционал ©, с учетом (2), получаем 
уравнения для определения функционала 21. Далее удобно ввести обозначение а = 1, в 
итоге получаем: 


де 2028 де 
-+6—>+(6, -В)0(х-х)-> -— 
&(Р) — В 22 (В, - В )6(х-ж) >. 
-2(/- 1.)6(=-х.)=0%И" (>). (15) 
В конечный момент времени и на обеих границах потребуем: 
вы 03 (16) 
08 08 
2 2. а 17 
дх х=хо дх х=хо о 


При этом вариация функционала принимает вид: 


(18) 


1х2 
де д 
Я] = —® — обийак = (УЛ, д, * . 
Гога В 
хо 


где градиент 


О а (19) 
Ох дх 


где х»+е|хо,х> |. Конкретное значение х»„ должно определяться в результате 
0,72 


анализа управляемости, а точнее — идентифицируемости [2]. 
Нетривиальное решение сопряженной задачи (15) — (17) возможно только при 
хе <х» <х|. Это объясняется наличием б-функции в уравнении (15) при линеаризо- 


ванном целевом функционале. Именно этот свободный член является источником 
нетривиального решения. При интегрировании уравнения (15) б-функция преобра- 


зуется к 9-функции отличной от нуля справа от х,. Следовательно, параболическая 


система (15) — (17) при приблизительном численном решении может иметь гаранти- 
рованное существенно ненулевое (нетривиальное) решение = справа от х.. 


Таким образом, функция (г) идентифицируема на 5 по целевому функцио- 
налу (6), если 5 = (1,1%) хх», хо <» <Х|. 
На основе градиента (19) организуется направленная итерационная коррекция 
функции А, (6) [2], [3]: 
ВАО = ВК -ЬокУЛ “О. (20) 
Здесь глубина спуска на каждой итерации вдоль выбранного направления ми- 
нимизации Бору. ^ (0) определяется числом Б^ по методу: 


если Л! </^, тогда БИ =ЫБ^, ы>Е 


тогда повторяется предыдущая 21 
с и (21) 
итерация, при В” =>Ь`; 0<Ё» <1. 
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Функция 4©(г)>0. Если @(Г) =1, то алгоритм (20) принимает вид алгоритма наи- 
скорейшего спуска [2], [4]. В выражении (21) принималось В =1.2, 55 =0.5. Функция 
@(г) регулирует направление спуска и определяется из условия не более 15% первого 
изменения идентифицируемой функции р, (Г): 


0.15. ° (Е) 


= (р) 


| (22) 


На рис. 1 приведены результаты идентификации коэффициента Д,(!) за 21 ите- 


рацию для реального процесса затвердевания стального расплава в изложнице на 
основе экспериментальных данных [7]. 


2. 


35 


25 


15 
0 4 8 +, с 


Рисунок 1 — Идентифицированный коэффициент 2, (г) 
Максимальное расхождение функции состояния составило всего лишь 
тах — 7, = 6°С. 


Это высокая точность моделирования процессов, описываемых математическими 
моделями с дифференциальными уравнениями параболического типа. 
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